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Abstract. We present direct measurements of the limb-darkened intensity profiles of the late-type giant stars 
HR5299, HR7635, and HR8621 obtained with the Navy Prototype Optical Interferometer (NPOI) at the Lowell 
Observatory. A triangle of baselines with lengths of 18.9 m, 22.2 m, and 37.5 m was used. We utilized squared 
visibility amplitudes beyond the first minimum, as well as triple amplitudes and phases in up to 10 spectral 
channels covering a wavelength range of ~ 650 nm to ~ 850 nm. We find that our data can best be described 
by featureless symmetric limb-darkened disk models while uniform disk and fully darkened disk models can 
be rejected. We derive high-precision angular limb-darkened diameters for the three stars of 7.44 mas ± 0.11 mas, 
6.18 mas ± 0.07 mas, and 6.94 mas ± 0.12 mas, respectively. Using the HIPPARCOS parallaxes, we determine linear 
limb-darkened radii of 114 Rq ± 13 R©, 56 Rq ±4R0, and 98 R© ±9R0, respectively. We compare our data to a 
grid of Kurucz stellar model atmospheres, with them derive the effective temperatures and surface gravities without 
additional information, and find agreement with independent estimates derived from empirical calibrations and 
bolometric fluxes. This confirms the consistency of model predictions and direct observations of the limb-darkening 
effect. 
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1. Introduction 

For a detailed understanding of stellar atmospheres and 
stellar evolution, it is of essential importance to obtain ac- 
curate observational estimates of stellar surface structure 
parameters on all scales. These parameters include diam- 
eters, limb-darkening profiles, photospheric asymmetries, 
and special features like hot spots. Model atmospheres 
are mainly constrained by observations of stellar spectra. 
Direct measurements of the limb-darkening profile can, in 
principle, provide an independent estimate of the temper- 
ature change with continuum opacity and, thus, an in- 
dependent observational verification. So far, the detailed 
intensity profile and the whole variety of additional sur- 
face structure parameters could only be observed in the 
case of the Sun. Today's interest, however, includes other 
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phases of stellar structure and evolution, fo r insta nce stars 
in late e volut ionary phases (Man duca e t al. 1977 , Sc holz fc 
Takeda [l98^ , Hofmann & Scholz |1998| , Jacob et al. ^OOOl) . 
These stars might exhibit asymmetric (e.g. Wittkowski et 
al. 1998| ) and even highly fragmented (e.g. Weigelt et al. 
1998 ) mass-loss events, which are believed to be triggered 
by the conditions on the stellar surfaces. Unfortunately, 
direct measurements of surface structure parameters are 
rare. While diameters have so far been obtained for sev- 
eral hundred stars with interferometric and lunar occulta- 
tion techniques, the second-order effect, limb-darkening, 
has been directly observed for only a very limited number 
of stars (Hanbury Brown et al 
Quirrenbach et al 



1974; Haniff et al. 1995 



1996; Hajian et al. 1998). Additional 



surface features have been detected on the apparently 
largest supergiants a Ori, a Sco, and a Her contributing 
with 5% to 20% to the total source fluxes at optical wave- 



Young et al. 2000) 



lengths (B uscher et al. |1990| ; Wilson ct al. |1992| ; Gil liland 
& Dupree 1996; Burns et al 



1997; TuthiU et al. 1997 
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The NPOI, located near Flagstaff, Arizona, is espe- 
cially designed for imaging of stars and their environments 
and is described in detail by Armstrong et al. ( 19981 ). 
The methods of "baseline bootstrapping" and "wave- 
length bootstrapping" (see Roddier 1988, Quirrenbach et 
al. 1996, Hajian et al. 1998) can be used in order to de- 
tect weak fringe contrasts, i.e. low visibility values, on re- 
solved stars by detecting the higher-contrast fringes on a 
chain of shorter effective spacings which comprise the long 
baselines. Hajian et al. (1998) demonstrated that boot- 
strapping with the NPOI astrometric subarray enabled 
the measurement of visibility values of a resolved stellar 
disk beyond the first minimum. By analyzing NPOI triple 
visibility products, they verified that the intensity profiles 
of a Ari and a Gas deviate from uniform disks due to the 
effect of limb-darkening. 

Here, we use NPOI's bootstrapping ability and apply 
an improved bias correction in order to utilize squared 
visibility amplitudes, in addition to triple amplitudes and 
closure phases, of three much fainter but well resolved late- 
type stars for spatial frequencies on both sides of the first 
minimum. Since there are three squared visibility ampli- 
tudes for each triple amplitude, more information can be 
used for the analysis of stellar intensity profiles. In an at- 
tempt to check the consistency of direct observations and 
model predictions of the limb-darkening effect as a func- 
tion of continuum wavelength, effective temperature, and 
surface gravity, we compare our multi-wavelength inter- 
ferometric data to a grid of Kurucz stellar model atmo- 
spheres. In addition, we determine high-precision limb- 
darkened diameters. 

2. Observations 

Interferometric observations of the three late-type giants, 
HR5299 (M4), HR7635 (K5), and HR8621 (M4), were 
performed with the NP OI us ing the configuration de- 
scribed by Benson et al. ( 1997 ). 

The center, east, and west siderostats of the astromet- 
ric subarray were used with effective apertures of 12.5 cm. 
They provide baselines with lengths of 18.9 m, 22.2 m, 
and 37.5 m at azimuths (measured east from north) of 
-67?5, 63?6, and 86?0, respectively. The three afocal 
beams of light are reduced in diameter and sent into vac- 
uum delay-lines for compensation of optical path differ- 
ences (OPDs) before they are combined pairwise using 
beamsplitters. Three of the afocal output beams corre- 
sponding to the three baselines are dispersed by a prism, 
focused by a lenslet array onto 32 optical fibers, and de- 
tected by avalanche photodiodes ( APDs) , covering a spec- 
tral range from 450 nm to 850 nm. The fringe packet is 
detected through modulation of the OPDs and is kept 
centered close to zero residual delay using the metho d of 
"group delay fringe tracking" (Armstrong et al. 1998 ). 

Table |l| lists names and characteristics of the observed 
stars together with the observing dates, the number of 
obtained scans, and the names and estimated diameters 
of the calibrator stars. During a scan of 90s, the pho- 



ton count rate for every channel is determined in eight 
temporal bins (synchronous with the delay line modula- 
tion), which sample a fringe every 2 ms. After each scan, 
a background measurement was taken on blank sky near 
the star. Immediately before or after each scan of a pro- 
gram star, a scan on one of the calibrator stars as specified 
in Table ^ was recorded. The calibrator stars were chosen 
to be located near the appropriate program stars on the 
sky. Their diameters were estimated using a calibration 



obtained by Mozurkewich et al. (1991) based on the ap- 
parent visual magnitude and the (i? — /) color index and 
are small enough so that possible errors in this estimate 
do not noticeably affect the calibration of our much larger 
resolved program stars. 

In order to compensate for detection noise bias terms 
(see Sect.^, incoherent (i.e. fringeless) data on several 
stars covering a range of apparent visual magnitudes were 
recorded on July 22, 2000, by moving the delay lines off 
the fringe packet. 

The signal-to-noise ratio of the measured visibilities 
decreases for spectral channels with shorter wavelengths 
owing to their narrower bandwidths (the channels are 
equally spaced in wavenumber), poorer seeing, and the 
red color of the observed stars. Therefore, only the 10 red- 
dest channels were used for the data analysis of HR7635, 
and the 5 reddest channels for that of the fainter stars 
HR5299 and HR8621. The central wavelengths of the 
10 reddest spectral channels are known to within about 
1% to be (852, 822, 794, 769, 745, 723, 702, 683, 665, 
649) nm, while their bandwidths range from '^31nm for 
the 852 nm channel to ~16nm for the 649 nm channel. 
These passbands mainly provide continuum observations 
and are not dominated by spectral features. The absorp- 
tion band/continuum ratio is relatively small for the K5 
star and the strong TiO bands at 671 nm and 714 nm are 
not covered by the spectral channels used for the M4 stars. 
However, for the M4 giants, some spectral channels are af- 
fected by molecular absorption bands, which will be taken 
into account. 

Figure |l| shows the obtained coverages of the uu-plane 
for HR5299, HR7635, and HR8621 based on all observa- 
tion dates and those spectral channels used for the data 
analysis. These spatial frequencies range up to a radius 
of ~300 cycles/arcsecond, which corresponds to a spatial 
resolution of 3.3 mas. 

3. Data reduction and calibration 

The raw data were processed, reduced, and calibrated as 



described in detail by Hummel et al. (1998). This process 
includes (1) calculation of the real iXi) and imaginary {Yi) 
parts of the complex visibility for each baseline i and each 
spectral channel by Fourier transform of the bin counts 
as a function of time, (2) calculation of the squared visi- 
bility amplitudes (|X^P) and amphtude (IV123I) and phase 
(0123, "closure phase") of the complex triple product, av- 
eraged over the 1 s intervals, (3) editing of the 1 s data, (4) 
further averaging over a scan of 90s, and (5) calibration 
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Table 1. Names (bright star catalog number HR, FK5 catalog number, common name) of the program stars, together 
with their properties (spectral type, apparent visual magnitude my, HIPPARCOS parallax TTtrig, bolometric flux i^Boi) 
and observational parameters (observing date, number of obtained scans, names of calibrator stars and their estimated 
diameters 0cai.)- Spectral Type, my, and TTtrig are taken from the HIPPARCOS catalogue (Perryman & ESA 1997), 
references for i^Boi are cited below. Note, HR7635 was classified as K5-M0 by Morgan & Keenan (1973), as K5 before. 



HR 


FK5 


Name 


Spectral 


mv 


7I"trig 


^Bol 


Ref. 


Observing 


#of 


Calibr. 


ecai. 








Type 


[mag] 


[mas] 


[10-1^ W/m^] 




Dates 


scans 


stars 


[mas] 


5299 


1368 


BY Boo 


M4.2III 


5.13 


7.01±0.66 


2.523±0.40 


a 


2000-07-07 


2 


FK5 527 


0.5 


















2000-07-13 


4 






7635 


752 


7Sge 


K5III 


3.51 


11.90±0.71 


2.792±0.14 


b 


2000-07-21 


6 


FK5 768 


0.3 


8621 




V416 Lac 


M4III 


5.11 


7.59±0.57 






2000-07-07 


2 


FK5 891 


0.3 


















2000-07-12 


2 


HR8494 


0.8 


















2000-07-13 
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References for bolometric fluxes: (a) Tsuji (1981), (b) Alonso et al. (199£) 
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Fig. 1. Plot of the obtained coverages of the uu-plane of the HR5299, HR7635 and HR8621 observations based on 
all observation dates (see Table |^) and on those spectral channels which were used for the data analysis, i.e. the 10 
reddest channels for HR7635 and the 5 reddest channels for the fainter stars HR5299 and HR8621. 



of the program star's squared visibility amplitudes, triple 
amplitudes, and closure phases by normalization with the 
corresponding smoothed (time kernel of 20 min) quantities 
of the calibrator stars. 

In step (2) \V\'^ is calculated using the unbiased esti- 
mator 

n 2 



=4 



7r/n 



sin(7r/n) 



< TV >2 ' 



(1) 



where n = 8 is the number of bins, N the total photon 
count rate in 2 ms, and the variance of the intensity 
caused by photon and detection noise. For stellar observa- 
tions, \V\'^ is compensated for background intensity by an 
additional factor < N >^ / < N - D where D is the 
background rate. The use of a squared quantity requires 
attention to the bias term (t| which is discussed in detail 
below. The triple product in step (2) is calculated using 
the unbiased estimator 

|Fi23|e*^-^ =< (Xi + iYi){X2 + iY2){X-, + zYg) > ■ (2) 

The phase (j)i23 of the complex triple product, formed by 
a triangle of baselines, is not corrupted by the phase noise 



caused by atmospheric turbulence (Jennison 1958). No 
bias correction is required for the triple product, since 
the noise from the three detector arrays receiving the sig- 
nal from each baseline is uncorrelated (e.g. Hummel et al. 
1998). In a triple with two short baselines and one long 
baseline resolving a stellar disk, as in the configuration 
used here, the triple amplitude has a higher signal-to- 
noise ratio than the squared visibility amplitude on the 
long baseline. This effect can be seen in our data shown 
below. It can be understood due to the fact that the low 
visibility on the long baseline is not squared but multiplied 
with higher amplitudes from the other two baselines in the 
triple. This effect was also confirmed by simulations of vis- 
ibility data based on Poisson noise. For the long baseline 
a visibility amplitude of 0.1 was assumed and for the two 
short baselines a visibility amplitude of 0.6. The signal-to- 
noise ratio of the squared visibility on the long (east-west) 
baseline was found to be 5 while that of the triple product 
was 20. 

Formal errors for the squared visibility amplitudes, 
triple amplitudes and closure phases were calculated based 
on the scatter of the 1 s samples. Calibration errors of the 
squared visibility amplitudes were estimated to be 7% for 
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HR7635 and 10% for HR5299 and HR8621 based on com- 
parisons of different scans and on calibrations with other 
cahbrator stars at shghtly farther distances in time and 
position. The total (formal and calibration) errors for the 
triple amplitudes and closure phases of HR 7635, HR 5299, 
and HR8621 were estimated to be 2, 1.5, and 2 times the 
formal errors, respectively. 



Noise bias compensation The use of the squared quantity 
\V\'^ requires a ttenti on to bias corrections as sho wn by 
Hummel et al. (|l998|) for NPOI data, by Colavita ( |l999|) 
for the case of t he Pa lomar Testbed Interferometer, and 
by Davis et al. (1999) for the Sydney University Stellar 
Interferometer. The noise bias term in Eq. |l| is equal 
to N in case of Poisson statistics (e.g. Shao et al. |l988| ). 
Since the NPOI detectors exhibit non-Poisson noise, due 
to after-pulsing of the APD s, the noise bias is estimated 
by erf = (Hummel et al. 1998 ). Here, Z"^ is the fringe 
amplitude floor estimated at four times the modulation 
frequency {X and Y are calculated at the temporal fre- 
quency k=l to select the component corresponding to the 
modulation frequency of the delay lines). However, de- 
spite the Z'^ bias compensation, positive visibility am- 
plitudes are observed for fringeless data and this bias, 
B := |yp (fringeless data, aj — Z^ , D = , was found to 
be a function of N and to be different for each spectral 
channel c and each spectrometer i. It was modeled with a 
power-law 



B{i,c) = ao(i,c)iV"i(*''=) 



(3) 



where parameters ao and ai were fitted to the fringeless 
data. This additional bias was compensated by using 



Z2 



(4) 



in Eq. 0. Obtained ranges for parameters oq and oi are 
[0.25,3.4] and [-1.48,-0.93], respectively. Thus, B is ap- 
proximately inversely proportional to N, i.e. the star's 
brightness. Consequently, its magnitude relative to the 
squared visibility amplitude is largest for low visibility val- 
ues of faint stars. Therefore, the compensation of this addi- 
tional noise bias is essential to obtain the very low squared 
visibility amplitudes of our faint program stars around and 
beyond the first minimum with a precision that allows an 
analysis of the limb-darkened intensity profiles. For in- 
stance, for our stars HR5299, HR7635, and HR8621, the 
total biases crj f/N^ for the bluest used channel and the 
37.5 m baseline amount to about 0.056, 0.063, and 0.063, 
respectively. Residuals of -0.004, -0.004, and ~0.003 re- 
main after the Z^ correction and are compensated by 
B. For comparison, the squared visibility amplitude of a 
fully darkened disk at the second maximum has a value of 
0.0074. 

The triple amplitudes and closure phases are not af- 
fected by this bias, as mentioned above. 



4. Data analysis and results 

Figures ||, |^ and ^ show all obtained squared visibility 
amplitude, triple amplitude, and closure phase data of 
HR5299, HR7635, and HR8621, respectively. Also shown 
are best fitting uniform disk models, fully darkened disk 
models, and Kurucz stellar model atmospheres as de- 
scribed below. 

The measured squared visibility amplitude of all three 
program stars as a function of increasing uw-radius de- 
creases monotonously towards a minimum, beyond which 
they increase. The low values are difficult to measure since 
they correspond to vanishing fringe contrasts. The mea- 
surement of almost zero squared visibility amplitude val- 
ues with acceptable error bars at the minima confirms the 
feasibility of the bootstrapping technique as well as the 
correctness of the photon and detection noise bias com- 
pensation. The observed functional form of the squared 
visibility amplitude is expected for a disklike object inten- 
sity distribution. The consistency of visibility values with 
baselines of different orientations excludes large deviations 
from circular symmetry. The absence of systematic varia- 
tions with smaller spatial frequencies excludes additional 
large-scale structures like circumstellar material with con- 
siderable intensity. 

The triple amplitudes show a similar behaviour for 
decreasing wavelength, i.e. increasing spatial resolution, 
with a minimum at a wavelength where the closure phases 
clearly exhibit a flip from to tt. The absence of intermedi- 
ate closure phases indicates object intensity distributions 
which are symmetric through reflection. Thus, for spectral 
channels with a bandwidth covering the location of the 
phase flip, intermediate values can occur. For HR8621, 
all recorded triple amplitudes and closure phases were be- 
yond the minimum and the flip, due to the star's large 
diameter in relation to the effective baseline lengths. The 
triple amplitudes and closure phases show a higher signal- 
to-noise ratio than the squared visibility amplitudes on 
the long baseline, as discussed in Sect. |[ 

The complex visibility of an astronomical object is re- 
lated to the object intensity distribution through a Fourier 
transform. Consequently, the object intensity distribution 
can in principle be directly reconstructed from interfer- 
ometric data using imaging techniques which effectively 
interpolate the limited coverage of the uw-plane. This was 
performed with NPOI data e.g. by Benson et al. (1997) for 
the double star Mizar A and by Hummel et al. (199?) for 
Matar. However, in order to obtain accurate estimates of 
physical parameters, model fits are a better choice. This 
applies especially to stellar disks since their imaging would 
require more resolution elements. Here, our data provide, 
roughly, two resolution elements across the diameter of 
the stars. The strength of the limb-darkening is related to 
the height of the second maximum of the visibility func- 
tion. The diameter is then determined by the locations of 
the minima of the visibility function and of the flip of the 
closure phases. 
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Fig. 2. Squared visibility ampli- 
tudes (left), triple amplitudes (mid- 
dle), and closure phases (right) of 
HR5299 and of best fitting models, 
for all obtained scans. The squared 
visibility amplitudes are plotted as 
a function of the uw-radius, and 
the triple amplitudes and closure 
phases as a function of the spec- 
tral channels' central wavelengths. 
The inset plots (left) provide an en- 
larged view of the low squared vis- 
ibility amplitude values around and 
beyond the first minimum. The x- 
symbols with error bars indicate the 
observations and their errors, the 
solid line the model values based 
on a uniform disk, the dashed line 
those based on a fully darkened disk, 
and the squares those based on the 
best fitting Kurucz stellar model at- 
mosphere (Kurucz 1993] ). The error 
bars include the formal errors as well 
as the calibration errors as described 
in Sect. 0. 
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Fig. 3. As Fig. I, but for HR7635. 
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Fig. 4. As Fig. I, but for HR8621. 
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Table 2. Best fitting diameters of HR 5299, HR7635, and 
HR8621 based on the models of a uniform disk (UD), 
a fully darkened disk (FDD) and the best fitting Kurucz 
model atmosphere (best K.) together with the obtained xt 
values. For HR5299, HR7635, and HR8621 the numbers 
of degrees of freedom are 150, 300, and 175, respectively. 
The quantity efinai denotes the final error of the diame- 
ter based on the best Kurucz model atmosphere due to 
the formal error (eformai) and calibration error (ecaiibr.) of 
the underlying data as well as on an error resulting from 
the choice of T^ft and logg (e„iodei)- The last four rows 
give the mean values and standard deviations of the UD 
and FDD diameters determined for each spectral chan- 
nel separately and corrected with factors derived from the 
expected model atmosphere. 





HR5299 


HR7635 


HR8621 




e 


xl 


e 


xl 


e 


xl 




[mas] 




[mas] 




[mas] 




UD 


6.82 


3.15 


5.67 


4.36 


6.25 


4.23 


FDD 


7.85 


1.38 


6.61 


2.26 


7.39 


1.67 


best K. 


7.44 


1.15 


6.18 


1.17 


6.94 


1.31 


*^formal 


0.03 




0.01 




0.03 




^calibr. 


0.10 




0.07 




0.10 




^modcl 


0.02 




0.02 




0.05 




Eflnal 


0.11 




0.07 




0.12 




corr. duo 


7.45 




6.11 




6.92 




a-{dvD) 


0.07 




0.02 




0.08 




corr. dFDD 


7.50 




6.13 




6.95 






0.06 




0.02 




0.08 





Uniform disk and fully darkened disk models Since the data 
as described above suggest a circularly symmetric disklike 
object intensity distribution /, models of a uniform disk 
(UD, J=lforO</i<l, / = otherwise, /i — cos 9 
being the cosine of the angle between the line of sight and 
the normal of the surface element of the star) and a fully 
darkened disk (FDD, I — fj) are used as a first approach 
to describe the data. In these cases, the visibility functions 
are given by 



Ui 



UD 



2 Ji(a^UD) 



^DD 



'^3/2(a^FDD) 
V ^ •''FDD 



(5) 



with a;uD,FDD — 7r0uD,FDD VM^ + i'^ a dimensionless 
spatial frequency (u, v: spatial frequencies [cycles/arcsec] 
as in Figures |l], ||, ^, 6: angular diameter of the star) 
and Ji and J3/2 the Bessel func tion o f first kind and orders 
1 and 3/2 (see e.g. Hestroffer |l997D . With V calculated 
for each baseline i € {1, 2, 3}, the squared visibility ampli- 
tudes |Vi|^, the amplitude of the triple amplitude IV1V2V3I 
and the closure phase 0i + (/>2 — 03 (0i = where Vi > 
and (f>i = TT where Vi < 0) can be derived. 

A minimization algorithm (simplex method) was 
applied in order to find the best fitting angular diameters 
GuD and 6fdd, using all available data, i.e. the squared 
visibility amplitudes, the triple amplitudes and the clo- 
sure phases. The derived diameters are shown in the first 



two rows of Table ||, together with the corresponding re- 
duced xl values. For HR5299, HR7635, and HR8621 the 
numbers of degrees of freedom are 150, 300, and 175, re- 
spectively. These xl values show that the UD model can 
be rejected and that the FDD model is a better descrip- 
tion of our data. The model functions are indicated by 
sofid lines (UD) and dashed fines (FDD) in Figs. |-|. The 
minima of the visibility and triple amplitude functions are 
well defined by the data. Before and in particular beyond 
the minima, our measured triple amplitudes for all three 
stars are significantly lower than the UD model values 
and slightly higher than the FDD model values. This in- 
dicates limb-darkened disks, less extreme than fully dark- 
ened disks, as predicted for late-type giants. Therefore, 
the specific limb-darkened profiles were investigated as de- 
scribed in the following paragraphs. 

Uniform disk and fully darkened disk models corrected for 
the effect of limb-darkening Before we present a direct fit 
of all our data to theoretical model atmospheres, we take 
an alternative approach that is based on results which are 
independent of a particular model atmosphere. 

Model atmospheres predict a decrease of the degree 
of limb-darkening with increasing wavelength, i.e. a tran- 
sition from near-FDD to near-UD continuum shapes, 
while the "true" limb-darkened diameter is wavelength- 

1977| , Hofmann & Scholz 



independent (Manduca et al. 



19981) . Since our data provide information at different 
spectral channels, they can be used to test stellar model 
atmospheres on their predicted wavelength dependence of 
UD and FDD diameters. 

Recently, Hestroffer ( |1997D discussed a limb- darkening 
law I\ = /z", with a a positive real number, as empir- 
ical brightness distribution function. This representation 
of the center-to-limb variation which uses only one limb- 
darkening coefficient includes the UD and FDD models 
and is very well suited to describe a wide range of dif- 
ferent realistic limb-darkening shapes (Hofmann & Scholz 
1998 ). Thus, in principle our data could be used to simul- 



taneously determine both the apparent limb-darkened di- 
ameter of our program stars and their limb-darkening pa- 
rameter a\ for each spectral channel. However, by means 
of Monte-Carlo simulations based on the actual coverage 
of the uu-plane of these observations and the claimed pre- 
cision, it was found that this determination is ambiguous. 

The employed method was first used by Hanbury 
Brown (1974), then applied by Quirrenbach et al. ( 1996[ ) 
and recently used for theoretical studies by Hofmann & 



Scholz ( |1998| ) and Davis et al. ( pOOOj ). The diameters based 
on UD and FDD models were derived for each spectral 
channel separately using only the squared visibility ampli- 
tudes up to the first minimum, since these data can well 
be described by UD and FDD models. The triple ampli- 
tudes and closure phases were not used since most of them 
contain visibility values beyond the first minimum, which 
do not fit UD and FDD models. The resulting UD and 
FDD diameters were multiplied by limb-darkening cor- 
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Fig. 5. Diameters based on uniform disk (UD, diamonds) and fully darkened disk (FDD, triangles) models for each 
spectral channel together with the Kurucz model atmosphere predictions (connected by solid lines). In addition, the 
limb-darkened (LD) diameters are shown which were derived by applying correction factors to the UD and FDD 
diameters. The mean LD diameters as quoted in Table ^ are indicated by the solid lines. The error bars represent 
the standard deviation of the obtained LD diameters. For comparison, the dotted lines indicate the limb-darkened 
diameters which were derived by the direct fit to a grid of Kurucz model atmospheres. 



rection factors to obtain the limb-darkened diameter of 
the star. Following the authors mentioned above, the cor- 
rection factors were derived as the ratios a;LD.o.3/2;uD,o.3 
and a;LD,o.3/a;FDD,o.3 with a;LD,o.3, a:^UD,o.3 = 2.0818, and 
2:fdd,o.3 = 2.3451 being the spatial frequencies where the 
squared visibility amplitudes |VldP, IHjdP and |T^fdd| 
(see Eq. H and||) equal 0.3. For the determination of xld, 
the Kurucz model atmosphere as described in detail in the 
following section was used with values for T^g and logy 
according to the spectral type of the star. 

The obtained diameters are plotted in Fig. ^ as a 
function of wavelength. The mean values and standard 
deviations of the corrected UD and FDD diameters are 
shown in the last four rows of Table ||. Our observed 
wavelength dependence of the UD and FDD diameters 
corresponds well with the model predictions. The small 
predicted deviations from a monotonous wavelength de- 
pendence are caused by effects of molecular absorption 
bands. Especially in the case of the brightest of our pro- 
gram stars, HR7635, our obtained diameters match the 
model predictions very well. Here, a larger number of spec- 
tral channels could be used and the effective temperature 
causes the transition from the UD to FDD model to oc- 
cur in the observed wavelength range. For all our program 
stars, no systematic deviations between observations and 
model predictions occur; which indicates the correctness 
of the atmosphere models used. Consequently, the derived 
limb-darkened diameters are, as required, independent of 
wavelength. The results based on the UD fit and those 
based on the FDD fit are consistent. 

The uncorrected UD and FDD diameters can be used 
for future comparisons with other model atmospheres and, 
furthermore, for future analyses taking additional infor- 
mation at other wavelengths into account. However, this 
approach cannot make use of our data at long base- 
lines and the corresponding triple amplitudes and closure 
phases. In order to provide more accurate limb-darkened 



Table 3. Comparison of minimum xt values obtained 
by fitting our HR5299 and HR7635 data to a grid of 
Kurucz model atmospheres with solar chemical abun- 
dances (Kurucz 1993 ) based on effective temperatures 
ranging from 3500 K to 4500 K and logg ranging from 
0.0 to 2.5. For each program star, the best xt values are 
marked by a box. For HR 5299 and HR 7635 the numbers 
of degrees of freedom are 150 and 300, respectively. 



HR 


logs/ 


0.0 


0.5 


1.0 


1.5 


2.0 


2.5 




















[K] 














5299 


3500 


1.25 


1.19 


1.15 


1.17 


1.19 


1.23 




3750 


1.19 


1.19 


1.19 


1.18 


1.18 


1.17 




4000 


1.21 


1.21 


1.21 


1.20 


1.20 


1.20 




4250 


1.27 


1.27 


1.26 


1.26 


1.26 


1.26 




4500 


1.42 


1.40 


1.38 


1.37 


1.36 


1.35 


7635 


3500 


1.50 


1.40 


1.35 


1.38 


1.47 


1.58 




3750 


1.29 


1.29 


1.29 


1.29 


1.27 


1.26 




4000 


1.18 


1.17 


1.17 


1.17 


1.18 


1.18 




4250 


1.20 


1.20 


1.19 


1.19 


1.19 


1.19 




4500 


1.26 


1.25 


1.24 


1.23 


1.23 


1.23 



diameters and to discriminate between different model as- 
sumptions a direct fit of all our data to Kurucz model 
atmospheres is described in the following paragraph. 

Comparison of our interferometric data with a grid of Kurucz 
stellar model atmospheres In an effort to compare different 
model atmosphere predictions based on a grid of effective 
temperatures and surface gravities with our interferomet- 
ric data, rather than assuming a particular model pro- 
file a priori, the best xt values based on different models 
were determined. In this way, theoretically predicted dif- 
ferences of the strength of the limb-darkening effect for 
different effective temperatures and surface gravities can 
be compared to direct measurements. 
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Kurucz (1993) tabulates monochromatic intensities 
/(I) and limb-darkening ratios for 17 val- 

ues of fi in 1221 frequency intervals ranging from 
9.09 nm to 160.0 /im, based on grids of model atmo- 
spheres for different chemical abundances. Here, his 
grid for solar chemical abundances and a microtur- 



bulent velocity of wtu 



2 km s was used (file 



"cfaku5.harvard.edu/grids/gridP00/ip00k2.pckl9"). The 
data is available for effective temperatures Toff ranging 
from 3500 K to 50000 K in steps of 250 K (for low Toff) 
and surface gravities \ogg (cgs) from to 5 in steps of 
0.5. For all our program stars, Tcff < 4000 K and logy < 2 
are predicted, based on their spectral types (see Table 
I). Thus, only models with Toff € [3500 K, 4500 K] and 
logy £ [0,2.5] were considered. 

NPOI passband-specific limb-darkening ratios 
Ic{fJ-) / Ic{l) were calculated by integrating the Kurucz 
data over each of the NPOI spectral channel's c € [1, 10] 
sensitivity functions. The NPOI spectral channels may be 
affected by molecular absorption bands, which is taken 
into account by the calculation of passband-specific limb- 
darkening model profiles. Compact photospheres were 
assumed, i.e. /c(m = 0) = 0. /c(l) was set to 1 since our 
measured visibility values are scaled to V{0) = 1 for each 
spectral channel separately. 

Following Davis et al. ( ^000 ), the model visibility val- 
ues were derived by numerical evaluation of the Hankel 
transform of the obtained tabulated intensity profiles 



/q /c(A^)Jo[a:LD(l - ^i'^y^^]fj.dfi 



(6) 



with XhD — Ti" ©LD Vv^ + v^ the dimensionless spatial fre- 
quency as described above, but now based on the star's 
limb-darkened angular diameter Old- No approximation 
of the tabulated model limb-darkened profiles by any limb- 
darkening law was used. The model squared visibility 
amplitudes, model triple amplitudes, and model closure 
phases were derived as described above. 

Based on each of the 30 considered Kurucz models (five 
values for and six values for logy), the wavelength- 
independent limb-darkened diameter, treated as the only 
free parameter, and the corresponding xt value were de- 
rived as described above using all our available data, i.e. 
the squared visibility amplitudes, the triple amplitudes, 
and the closure phases. 

For our three program stars, the lowest obtained xt 
values are listed in Table |^, together with the obtained 
best fitting limb-darkened diameters. Table || shows all 
resulting xt values for HR5299 and HR8621, where the 
lowest xt values are marked by a box. 

The occurrence of minimum xt values larger than 1.0 
might, in principle, be caused by optimum parameters ly- 
ing in between our grid points, a wrong model assumption, 
an underestimation of the calibration errors, or system- 
atic calibration errors. An underestimation of the calibra- 
tion errors leading to total errors underestimated by only 
7% and 8% is most likely to be the main cause for the 



deviations from unity in the cases of HR5299 and HR 
7635 since these errors can only be roughly estimated (see 
Sect. H). A considerable part of the larger deviation from 
unity, a value of 1.31, in the case of HR 8621 might also be 
caused by an effective temperature lower than 3500 K, i.e. 
an incorrect model assumption, or systematic effects dur- 
ing the calibration process. The stellar atmosphere model 
was adopted as the best fit to our data for HR 5299 and 
HR 7635. The xt values were analyzed as a function of 
Tcff and logy for these program stars only, as follows. 

For these program stars, significantly different xt val- 
ues are obtained for different model parameters, for exam- 
ple, in the case of HR5299 we obtain values between 1.15 
and 1.42. At higher temperatures, differences for varying 
Tcff are larger than for varying logy, because of the lat- 
ter's lesser effe ct on the limb-darkened profile (see e.g. 
Manduca et al. |1977D . 

To take the deviations of the xt values from unity 
into account, the values in Tab. ^ were normalized to 
unity at the minimum for the following analysis. Using 
this method, the total assumed errors in our data are in- 
creased by common mean factors of 1.07 and 1.08, ne- 
glecting that calibration errors depend on the value of the 
visibility. It was verified that due to the smallness of this 
correction the results obtained are still valid and that it 
is insignificant whether the total data errors are rescaled 
or just the calibration errors. 

Near the minimum, the function is expected to be 
a quadratic function of each of the varied parameters. 
Therefore, for each star, a parabola was fitted to the 
values as a function of Toff with fixed best-fitting logy 
and as a function of logy with fixed best-fitting Toff. Here, 
more digits were used than shown in Table |^. The most 
likely values for the parameters Toff and logy can be es- 
timated by the locations of the minima of the parabola. 
Assuming purely Gaussian noise, the corresponding 1 a 
errors can be estimated as the variation in the parame- 
ters which will increase the normalized total values by 
1 from its value at the minimum of the fitted parabola 
(see e.g. Bevington & Robinson 1991). For HR5299 and 
HR7635, the x^ values as a function of Toff as well as of 
logy match a parabola very well. This confirms that ef- 
fects due to systematic calibration errors or an incorrect 
model assumption are not of considerable size. However, 
small additional errors due to these effects cannot be ruled 
out and are not included in the error analysis presented 
here. For HR5299, the xt values in Table ^ as a function 
of Toff extend to the minimum but not beyond. However, 
the one-sided xt values fit a parabola with a minimum at 
3520 K very well, confirming that our grid point at 3500 K 
is in fact close to the minimum. The derived most likely 
values and the errors are shown in Table ^ and are com- 
pared to independent estimates in Sect. ^. 

For all our program stars, the best xt values derived 
here are significantly better than those based on the FDD 
model (see Table ||). This effect is most noticeable in the 
case of HR7635, due to the higher signal-to-noise ratio of 
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the data and the higher Tcff resulting in hmb-darkening 
that is not as close to the FDD case as for the cooler stars. 

The limb-darkened diameters corresponding to the 
best fitting models are shown in Table Errors egnai were 
derived based on the formal errors and calibration errors 
mentioned in Sect. ^ and on those due to the choice of 
Teff and logg. For the latter error, the standard devia- 
tions of diameters based on models with Teff and logg 
values within their error bars were taken. For HR8621 all 
30 considered models were included. 

The squared model visibility amplitudes, model triple 
amplitudes, and model closure phases obtained with the 
best fitting Kurucz model atmosphere are indicated by the 
squares in Figs. ^-|^. They coincide well with our measured 
data and describe them considerably better than the UD 
and FDD models. The observed data of HR8621 differ 
slightly from the model values which might be explained 
by calibration errors or a wrong model assumption as men- 
tioned above. The absence of further systematic deviations 
between the model predictions and our data confirms that 
extended photospheres with I{fi = 0) > need not to be 
considered and that the width of the NPOI spectral chan- 
nels does not noticeably affect our analysis. 

Deviations from featureless symmetric disks Although, as 
discussed above, our data indicate a featureless symmet- 
ric disk, a formal search for solutions with an asymmetric 
object intensity distribution was performed in order to 
validate this interpretation. All our data were fitted to 
a limb-darkened profile as derived before but with an el- 
liptical shape. Starting values for the position angle of 
the major axis were chosen from 0°to 180°in steps of 
10°and for the axis ratios (major axis/minor axis) two 
starting values of 1.05 and 1.1 were considered. For the 
stars HR5299, HR7635, and HR8621, slightly better xl 
values than for the symmetric intensity profile were de- 
rived with axes ratios of 1.02, 1.01, and 1.01, respectively. 
These small deviations from circular symmetry are not 
significant for an asymmetric object intensity distribution, 
but can be caused, for instance, by systematic calibration 
errors. However, this study does exclude elliptical object 
intensity distributions with larger axis ratios. 

In order to estimate whether additional hot spots, as 
they were found on the surfaces of a Ori, a Sco, and a Her 
(see Sec. |^) could be detected in our data, model squared 
visibility amplitudes, triple amplitudes, and closure phases 
were calculated for one example based on our best fitting 
limb-darkening model with one additional hot spot. The 
spot's intensity was chosen to be only 2% of the star's 
intensity, which is clearly less than that of the spots on 
a Ori, a Sco, and a Her. The spots were assumed to be 
unresolved and to have a separation of half of the star's 
limb-darkened radius. Two position angles were consid- 
ered, namely that of the preferred direction of the uw-plane 
coverage and that perpendicular to it. Figure |^ shows the 
model predictions, indicating that the existence of such 
a spot would significantly affect the closure phases and 



the triple and visibility amplitudes around their minima. 
The different triple and visibility amplitudes might be 
modeled by another stellar diameter and limb-darkening 
profile, but the occurrence of closure phases significantly 
different from values of and tt, however, could only be 
explained by an asymmetric intensity distribution. Our 
data is not consistent with such a noticeable asymmetry. 
Consequently, it can be concluded that the existence of an 
unresolved single spot on the surfaces of our program stars 
with an intensity at least as high as studied above is highly 
unlikely. A resolved spot is unlikely, too, since its intensity 
contribution would be higher. However, by an analysis of 
the closure phases we cannot rule out a centered spot. 

5. Discussion 

Based on the data analysis described in Sect. || it fol- 
lows that all our interferometric data are consistent with 
featureless circularly symmetric limb-darkened disks. The 
method of comparing all our interferometric data to a grid 
of Kurucz model atmospheres allows the determination of 
the effective temperature and surface gravity without ad- 
ditional information in the cases of HR 5299 and HR 7635 
and provides an accurate estimate of the limb-darkened 
diameter Old- 

Independent estimates for Toff and \ogg can be ob- 
tained by empirical calibrations of the spectral type. 
Additionally, Toff and the linear limb-darkened diameter 
D can be derived by using the bolometric flux -FboI (see 
Table |l|) and the HIPPARCOS parallax TTtrig, together 
with our value for Old- The errors are dominated by the 
uncertainties of i^Boi and TTtrig- As a result, variations of 
Old have little impact on Tog and this estimate of Toff 
can be regarded as sufficiently independent of our deter- 
mination by the direct fit to Kurucz model atmospheres 
as well. 

Table ^ lists the results obtained by the direct fit to 
Kurucz model atmospheres together with the independent 
estimates. The error estimates for T^s and log 5 are based 
on the analysis of Table ^ as described above. 

Figure]^ compares the two derived T^s values, obtained 
by the direct fit to a grid of Kurucz model atmospheres, 
and by Old and Tboi, with different empirical calibrations. 
For HR8621 the data quality is, as mentioned above, not 
sufficiently high to obtain Teff and log g by means of the di- 
rect fit to Kurucz model atmospheres. For HR5299, both 
the Toff value and log g value derived by this fit to Kurucz 
model atmospheres are well consistent with the indepen- 
dent estimates. For HR7635, the obtained Toff value is 
higher than that of the independent estimates and con- 
sistent only within ^ 2 a. The value obtained for logg is 
consistent with the empirical calibration. 

Our derived values for Old are generally consis- 
tent with earlier determinations of uniform disk diame- 
ters corrected for limb-darkening. They are available for 
HR5299 (7.0 mas ±0.3 mas by Dyck et al |l996i |l998h 



and for HR7635 (7.4 mas ±0.2 mas by Hut ter et al 



5.5 mas ± 0.5 mas by Dyck et al. |l996i |l998D . For HR 7635 
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Fig. 6. Impact of hot spots on model 
data for the example of HR7635, 
scan no. 5. The triangles denote the 
best fitting limb-darkening model 
with one additional hot spot. The 
spot's intensity relative to that of 
the star is 2%, its separation is 
^ld/4, its position angle is (top) 
along the preferred direction of the 
uv plane coverage (3°) and (bot- 
tom) perpendicular to it (93°). For 
comparison, the x-symbols with er- 
ror bars and the squares denote the 
measurements and the best fitting 
limb-darkening model, respectively, 
as in Fig. |3[ 



Table 4. Best fitting limb-darkened diameters Old, effective temperatures, and surface gravities as derived by the 
direct fit of our interferometric data to a grid of Kurucz model atmospheres as described above (Tables |^ and 
Using Old together with the HIPPARCOS parallax TTtrig and the bolometric flux FboI (see Table |l|) the linear stellar 
diameter D (Old, '"'trig) and the effective temperature Toff (Old, ^Boi) are derived. For further comparison, the effective 
temperatures T cff (Sp .T.) and surface gravities logg (Sp.T.) derived from empirical calibrations of the spectral type 
(Schmidt-Kaler |1982| ) are shown. 



HR 


Old 

[mas] 


Teff [K] 




D (0LD, ""trig) 

[R©] 


left (6lD, -Fboi) 

[K] 


Toft (Sp.T.) 
[K] 


\ogg (Sp.T.) 


5299 


7.44±0.11 


3520±190 


1.3±0.4 


228lf3 


3420±160 


3410 


<1.3 


7635 


6.18±0.07 


4160±100 


0.9±1.0 


112^8 


3850± 70 


3950 


1.7 


8621 


6.94±0.12 






19711? 




3430 


<1.3 



Alonso et al. 200C derived a limb-darkened diameter of 
6.12 mas ±0.2 mas by means of the infrared flux method 
(IRFM), wich is in very good agreement with our value. 

Our derived linear limb-darkened radii for the three 
stars of 114R0±13Ro, 56R0±4R0, and 98R0±9R0 
are in good agreement with those obtained with the em- 
pirical calibration for M giants by Dumm & Schild ( 1998 ) 
based on the HIPPARCOS parallaxes, V magnitudes, and 
V — I color indexes, which are 119 R0, 60 R0, and 97 R0, 
respectively. 

The circular symmetry of our observed object inten- 
sity distributions is expected because at optical wave- 
lengths only the surfaces of the stars themselves are 
observed rather than additional circumstellar envelopes 
where asymmetric morphologies were discovered. These 
asymmetries in the envelopes can be caused by, e.g., rota- 
tions so slow that they do not observably affect the star's 
shape. The absence of additional surface features as ob- 
served on the surfaces of the apparently largest super- 
giants might be explained by the higher surface gravities 
of our program giant stars. 



6. Summary 



Featureless symmetric limb-darkened stellar disks provide 
good fits to our NPOI interferometric data of IIR5299, 
HR7635, and HR8621. We are able to discriminate be- 
tween model atmospheres with different effective temper- 
atures and surface gravities. We find that our interfero- 
metric measurements and stellar model atmosphere pre- 
dictions by Kurucz ( 1993| ) of the limb-darkening effect are 
consistent. We obtain high- precision (l%-2% accuracy) 
limb-darkened angular disk diameters and derive linear 
radii and effective temperatures using the HIPPARCOS 
parallaxes and bolometric fluxes reported in the litera- 
ture. With reduced noise terms and the upcoming simul- 
taneous combination of six beams at NPOI we will be able 
to obtain even more precise limb-darkened diameters of a 
much larger number of stars. Furthermore, interferomet- 
ric data will soon allow the discrimination between model 
atmospheres with different effective temperatures, surface 
gravities and even chemical abundances with higher preci- 
sion than in this first attempt. Thus, further observational 
constraints for model atmospheres will become available, 
in addition to observations of stellar spectra. 
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values derived from the limb-darkened angular diameter 
0LD and the bolometric flux Tboi- 
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